Structural precursor to the metal-insulator transition in V2O3 
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The temperature dependence of the local structure of V2O3 in the vicinity of the metal to insulator 
transition (MIT) has been investigated using hard X-ray absorption spectroscopy. It is shown that 
the vanadium pair distance along the hexagonal c-axis changes abruptly at the MIT as expected. 
However, a continuous increase of the tilt of these pairs sets in already at higher temperatures and 
reaches its maximum value at the onset of the electronic and magnetic transition. These findings 
confirm recent theoretical results which claim that electron-lattice coupling is important for the 
MIT in V2O3. Our results suggest that interactions in the basal plane play a decisive role for the 
MIT and orbital degrees of freedom drive the MIT via changes in hybridization. 

PACS numbers: 71.30.+h, 61.10.Ht, 71.27. +a, 78.70.Dm 



Introduction 

The metal to insulator transition in V2O3 is a long 
standing yet not completely resolved issue. Originally 
being discussed as the classical example of a Mott- 
Hubbard transition^ this interpretation has to be ques- 
tioned due to the intimate interplay of structural, elec- 
tronic and magnetic changes in the phase diagram of this 
compound.- Therefore it is not clear whether electronic 
correlations can completely account for the MIT in V2O3 
or if structural effects play an important role. The rel- 
evance of possible electron-lattice interactions has been 
stressed by recent theoretical worki&i 
In an earlier paper— we have shown that a detailed under- 
standing of the interplay of structural and electronic de- 
grees of freedom requires both information on the global 
structure as measured by X-ray diffraction (XRD) and 
of the local structure of the compound. Interatomic dis- 
tances and symmetries can deviate locally from the global 
lattice symmetry, therefore indicating the presence of an 
order-disorder component of the transitions or of struc- 
tural fluctuations. Our earlier results demonstrate that 
the two insulating phases of V2O3, the paramagnetic in- 
sulating (PI) and the antiferromagnetic insulating (AFI) 
phase, are structurally very similar on a local scale and 
therefore indicate that there might exist a common route 
for both the transitions from the paramagnetic metallic 
(PM) phase to the PI and the AFI phase, respectively. 
In this letter, we investigate the evolution of the local 
structure of V2O3 in the vicinity of the PM to AFI tran- 
sition to shed light on the detailed structural changes 
associated with the MIT. The results of our tempera- 
ture dependent extended X-ray absorption fine structure 
(EXAFS) measurements show a rather complex behav- 
ior of the local structure, which differs along different 
crystallographic axes. In particular, we observe a struc- 
tural precursor to the MIT which consists of character- 
istic changes of nearest neighbor distances for the first 
time. 



Experimental 



EXAFS measurements were performed on a V2O3 single 
crystal grown by chemical transport using TeCLj as trans- 
port agent. The crystal was mounted on a thin glass plate 
and reduced in thickness by lapping and polishing the 
sample down to approximately 12 /im which corresponds 
to roughly two absorption lengths. This treatment al- 
lowed to detect the X-ray absorption in a transmission 
geometry, thereby avoiding damping of the signal due to 
selfabsorption effects which would have been present in 
fluorescence mode— The sample was oriented so that the 
angle between the polarization vector E of the incom- 
ing X-rays and the hexagonal c axis of the crystal (chex) 
could be changed by rotating the sample around the sur- 
face normal. Measurements were made with E parallel 
and perpendicular to c^ex ■ This procedure allows to sep- 
arate scattering paths along and perpendicular to the 
hexagonal c axis as detailed in ref. |5| 
The electrical resistivity was measured using a standard 
four-probe DC technique. For each measurement the cur- 
rent was reversed to eliminate offsets. Contacts were 
made by attaching platinum wires (diameter: 40 /nm) 
to the sample by silver epoxy. Metallic behavior was ob- 
served in the high temperature region, and a sharp jump 
in the resistance of more than three orders of magni- 
tude at a temperature of 140 K signaled the MIT. 7 The 
slightly reduced transition temperature indicates that a 
small vanadium deficiency occurred during the growth 
of the single crystal. The magnetic transition from PM 
to AFI occurred in a temperature region from 140 K 
down to 130 K as determined by magnetic susceptibility 
measurements using a Quantum Design superconducting 
quantum interference device (SQUID) magnetometer. 
EXAFS measurements at the vanadium K edge (5465 
eV) were performed at beam line X23B of the National 
Synchrotron Light Source at Brookhaven National Lab- 
oratory. Spectra at different temperatures above, at, 
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FIG. 1: (color online). EXAFS spectra of the V K edge of 
V2O3 in the vicinity of the MIT at T=140 K. The spectra 
in a geometry with E parallel to Chex are almost identical 
within one phase and display an abrupt change at the MIT 
(a) . Only smooth changes are observed in a setup where E is 
perpendicular to Ch ex (b). 



and below the transition temperature were taken using a 
closed cycle helium cryostat. The sample was held under 
vacuum to avoid water condensation and to reduce heat 
transmission. As no temperature sensor was available in 
the direct neighborhood of the sample, the resistivity of 
the crystal was monitored in situ during the EXAFS ex- 
periment in the geometry with E perpendicular to Chex- 
This allowed for a direct observation of the metal to insu- 
lator transition itself. The electronic transition was ob- 
served at a temperature of 130 K at the cold finger of the 
cryostat. The temperature difference of -10 K compared 
to the ex situ resistance measurements can be attributed 
to a temperature gradient between the cold finger and the 
sample. Temperatures given in the paper are corrected 
accordingly. 

The crystal has been examined by XRD measurements 
after the EXAFS experiments were finished. The inves- 
tigations, both in 6/ 20 and in Laue geometries, con- 
firmed that the MIT did not destroy the single crys- 
tal, as the diffraction patterns were those expected for 




FIG. 2: (color online) . Room temperature structure of V2O3. 
Small (large) spheres represent vanadium (oxygen) ions. Only 
those oxygen ions which form an octahedron around the cen- 
tral vanadium ion Vo are shown. Vanadium ions which have 
opposite spin in the AFI phase are drawn in different colors. 



a good single crystal and no signs of twinning could be 
observed. In addition, spectra above the transition were 
measured before and after the crystal passed through the 
MIT during the EXAFS experiment. No differences were 
observed in these spectra, ruling out the influence of pos- 
sible cracks on the measurements as well. Typcial spec- 
tra are shown in fig. ^ and display that a good signal to 
noise ratio could be achieved at least up to a wavevector 
of k = 14 A -1 (k being the wavevector of the outgoing 
photoelectron). 

Standard EXAFS analysis using the ATHENA^ and 
IFEFFIT 9 programs was performed to extract path 
lengths and other parameters from the spectra. The 
experimental spectra were normalized to the edge step 
and a smooth, atomic-like background was subtracted 
by minimizing the low-i? contribution in the Fourier- 
transformed dataiiS The EXAFS equation in cumulant 
expansion as defined in Ref. ^] and implemented in the 
IFEFFIT program was used up to the third cumulant C3 
to fit the experimental data to model spectra calculated 
with FEFF8.2ii2 For the Fourier transforms a window in 
k space from 2 A -1 to 13 A -1 was used and the spectra 
were weig hted by fc 2 5 . Fitting was carried out in real 
space in a range from 1.2 A to 4.1 A (E \\ Chex) an d 
from 1.2 A to 3.65 A (E J_ c/j ex ), respectively. The clus- 
ter of ions used for the calculation of the model spectra 
was derived from the lattice parameters and atomic posi- 
tions resulting from XRD measurements for the trigonal 
structure of V2O3 at room temperature published in Rcf. 
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TABLE I: Definition of the next neighbor V-V distances in 
the fitting model. The free parameters are dvl for the change 
of the Vo-Vi distance and dv2 for the change of the Vo-V2a 
distance, respectively. The changes AR of the remaining V- 
V distances used in the fits with with E \\ c^ex and with 
E _L Chex are set as multiples of these values to represent the 
values calculated from the trigonal and monoclinic structure 
as measured by XRD (right columns) when dvl and dv2 take 
on the respective values of dvl = 0.04 A and dv2 = 0.11 A. 



The structure refinement was performed by varying the 
interatomic distances to all vanadium ions up to a shell 
radius of 4.3 A, which corresponds to the V ions shown 
in fig. In addition, backscattering from the six oxy- 
gen ions that form the coordination octahedron of the 
absorber was accounted for. Inclusion of the second 
backscattering shell of oxygen ions improved the quality 
of the fit, but had no influence on the fitting parameters 
concerning the ions shown in Fig. [21 At the transition 
from the paramagnetic metal to the antiferromagnetic 
insulator, the lattice symmetry reduces from trigonal to 
monoclinic. Comparing the interatomic distances of the 
trigonal and the monoclinic structure as calculated from 
the data published in Refs. EH and lillTlill and sum- 
marized in the two rightmost columns of table [IJ it turns 
out that the transition increases the distance of the neigh- 
boring vanadium ions along the hexagonal c axis (Vo and 
Vi) by 0.04 A. In addition, the threefold symmetry in the 
plane perpendicular to Ch ex (the basal plane) is broken 
and one of the three vanadium neighbors in this plane 
(V2a) increases its distance to the absorber (Vo) by 0.11 
A, while the distance to the two other V2 ions remains 
nearly constant. 

To allow for this monoclinic distortion in the fitting 
model, the distance between the absorber and these V 
neighbors (Vi for the spectra with E || Ch ex and V2 a for 
the spectra with E _L Chex) was allowed to vary freely, 
while the distances to the remaining V ions were given 
as multiples of this value as defined in table |U This al- 
lowed to describe the changes of all V-V distances with 
a single parameter for each orientation (dvl or dv2). In 
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TABLE II: r-factors and reduced chi-squared (xl) for the fits 
presented in figs. 13151 

the same way the positions of the O ions were split up 
as multiples of dvl and dv2, respectively. An additional 
isotropic change of the V-0 distances was allowed to ac- 
count for a possible change in the oxygen positions due 
to the vanadium deficiency of the crystal. A damping 
term (Debye- Waller factor) was assigned to each neigh- 
boring shell, and for the nearest vanadium neighbors the 
third cumulant was used as an additional free parameter. 
Altogether, 12 free parameters were used to fit the 16.7 
independent points of each spectra. 



Results 

All measured spectra could be fitted successfully with 
a reliability factor of less than 0.02, indicating that the 
model was well suited to describe the data. For temper- 
atures slightly below and above the MIT, the Fourier- 
transforms of four spectra and the respective fits are 
shown in fig. |3J The values for the reliability factors and 
reduced chi-squared for the fits are compiled in table ITT1 
Figures 01 and |S] show the fit results for the scattering 
paths to the next vanadium neighbors along (Vi) and 
perpendicular (V2, V2 a ) to the hexagonal c axis. The ge- 
ometry used in the measurement allowed for a separation 
of these paths: Backscattering from the Vi neighbor only 
occurs in the spectra with E \\ Chex, while the V2 and 
V2a backscatterers only contribute to the spectra with 
E _L Chex- The Fourier transformed spectra show a well 
separated structure located between 2 A and 3 A (see 
fig. [SJ. This structure is attributed to single scattering 
contributions either of the Vi or the V2, V2 Q scatter- 
ing paths, only. Therefore no correlations between fit 
variables of these paths can occur, although the distance 
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FIG. 3: (color online). Typical examples of the Fourier trans- 
formed EXAFS spectra of V2O3 at temperatures close to the 
MIT in a geometry with E parallel to Chex (a) and with E per- 
pendicular to Chex (b). Solid and dash-dotted lines represent 
the amplitude and imaginary part of the Fourier transformed 
data, dashed and double-dotted dashed lines show the fitted 
model spectra. Data and fit coincide well in the fitting range 
(marked by vertical dotted lines) as can be seen from the dif- 
ference spectra of the imaginary parts of data and fit in the 
respective lower panels. For this figure, the scattering phase 
shift of the central V ion has been corrected, which basically 
leads to a shift of the fourier transformed spectra and the 
fitting ranges to the right by ~ 0.3 A. 



from the absorber Vo to all of these neighbors is similar. 
Figure Ufa) clearly shows that the distances between the 
next V neighbors along the hexagonal c axis correspond 
to the values expected from XRD results: The distance in 
the AFI phase at low temperatures is about 0.04 A larger 
than in the high temperature PM phase. The change is 
abrupt as a function of temperature and occurs in a tem- 
perature window corresponding to that of the change in 
magnetic susceptibility. The measurement uncertainty of 
the V0-V1 distances amounts to « 0.01 A. The Debye- 
Waller factor a 2 , shown in Fig. 02a), exhibits a sharp rise 
at the onset of the MIT, indicating either large fiuctua- 
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FIG. 4: (color online). Temperature dependence of the dis- 
tances of the V0-V1 (a) and Vo-V~2a (b) scattering paths as re- 
sulting from the fits to the EXAFS spectra (dots). Solid lines 
represent the distances expected from XRD. The dash-dotted 
line in (b) shows the result of the ex situ resistance measure- 
ment. The shaded temperature range marks the range of the 
magnetic transition. Dashed lines through the fit results pro- 
vide a guide to the eye. 



tions or a static distribution of widely varying distances 
while the magnetic transition occurs. 
The next neighbor distances in the basal plane display 
a very different behavior. First of all the distance from 
the absorber ion Vo to the V2 a neighbor is larger than 
expected from XRD in the high temperature PM phase 
[see Fig. E[b)], i.e. larger than the V0-V2 distances. This 
break of the threefold symmetry expected from the long 
range structure probed by XRD corresponds to a mono- 
clinic symmetry. The elongation of the Vo-V"2a distance 
, which could be determined with an accuracy of « 0.02 
A, measures the size of this monoclinic distortion. In the 
PM phase it is reduced to about 30 % of the distortion 
in the AFI phase, as was first observed in Ref. [ijj The 
transition to the full monoclinic distortion in the basal 
plane is much broader in temperature than that observed 
for the V-V distance along the hexagonal c axis. The on- 
set occurs in a region well above the electronic and mag- 
netic transitions. This is consistent with the fact that 
the Debye- Waller factor of the in-plane V neighbors [Fig. 
0V)] does not show a significant change at the MIT. 
Third cumulants (C3 or a^) which account for the 
anisotropy of the scattering paths to the next V neigh- 
bors where found to have a significant influence on the 
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FIG. 5: Temperature dependence of the Debye- Waller factors 
a 2 for the Vo-Vi (a) and Vo~V2,2a (b) scattering paths as re- 
sulting from the fits to the EXAFS spectra. Like in Fig. 0] the 
shaded temperature range marks the range of the magnetic 
transition. The dashed line through the fit results provides a 
guide to the eye. 



spectra in the low temperature AFI phase, but are neg- 
ligible in the PM phase. The third cumulants show only 
a weak temperature dependence in the AFI phase and 
adopt values of w —0.0003 A 3 for the Vo-Vi scatter- 
ing path and « -0.0007 A 3 for the V -V 2 and V -V 2a 
scattering paths, respectively. The error is estimated to 
approximately ± 25 %. 

Distances between the central vanadium ion Vo and the 
oxygen ions forming its coordination octahedron did not 
deviate from the values expected from XRD within error 
bars, i.e. these vanadium-oxygen distances remain nearly 
constant throughout the whole temperature range. 
Summarizing, we can characterize the structural evolu- 
tion at the metal to insulator transition of V2O3 as fol- 
lows: Structural changes along the hexagonal c axis are 
restricted to a sharply bounded temperature range given 
by the MIT. Perpendicular to Ch ex , however, the tilt of 
the c-axis V pairs, indicated by the elongation of the Vo- 
V 2a bond, starts to increase already far above the MIT. 



Discussion 



Precursor effects to the MIT in V2O3 have already been 
observed in electrical resistance under pressure^ and 



very recently in the sound velocity measured for thin 
films The results of our local structure determination 
show for the first time that structural changes preceed 
the changes in transport and magnetic properties. 
The analysis of the polarization dependent measurements 
reveals that two distinct movements with different tem- 
perature dependencies have to be distinguished: The po- 
sitioning of the vanadium ions within the basal plane 
and the change of the V-V pair distance perpendicular 
to the basal plane, i.e. along the hexagonal c-axis. The 
former movement determines the rotation of these vana- 
dium pairs, ranging from an alignment parallel to C} lex in 
the PM phase to a tilt of roughly 2° away from this axis 
in the AFI phased 

Following Tanaka^ hybridization between the a\ g and 
the orbitals increases on rotating the Vo-Vi pairs 
out of the direction of the c-axis. Depending on the hy- 
bridization strength, two minima of the free energy exist, 
which are characterized by two different equilibrium dis- 
tances between the two nearest V neighbors along Ch ex - 
While the electron-lattice coupling is weak when the V-V 
pair is aligned parallel to the hexagonal c-axis, increased 
hybridization on tilting the pair renders the orbital con- 
figuration almost energetically degenerate. The so mod- 
ified interactions favor a new ground state configuration 
with a reduced occupation of the a\ g orbitali22i2I In this 
configuration, the electron-lattice coupling is strong and 
an enlarged V pair distance is expected. 
We observe a continuously increasing tilt, starting well 
above the MIT and continuing even within the AFI 
phase. Along with it, the hybridization strength increases 
as just pointed out. At a certain hybridization strength, 
the system snaps into the second minimum of the free 
energy, which corresponds to the enlarged pair distance. 
Therefore it is plausible that the Vrj-Vi distance ex- 
pands primarily after the tilt of the pair has adopted 
its maximum value, what is exactly the behavior we ob- 
serve experimentally. Our measurements thereby sub- 
stantiate the scenario proposed by Tanaka which does 
not regard the metal- insulator transitions in V2O3 (both 
the PM to AFI and the PM to PI transitions) as usual 
Mott transitions^. The fact, that on a local scale the 
trigonal symmetry is not fully recovered in the metallic 
phase, suggests that the trigonal and monoclinic struc- 
ture are very close in energy, allowing fluctuations be- 
tween the two. We want to point out that the present 
results are consistent with our earlier findings^ that the 
trigonal symmetry can be broken on a local scale while 
the global symmetry remains trigonal. 
In step with the separation of the two neighboring V ions 
along the hexagonal c axis an "umbrella like" distortion 
of the oxygen coordination octahedron around each V 
ion occurs^ Thus, the trigonal distortion of the coordi- 
nation octahedron is directly determined by the distance 
between the next V neighbors along Ch ex ■ The increased 
distortion of the oxygen octahedron is responsible for the 
stabilization of the insulating stated since it results in a 
large shift of (dynamical) spectral weight. 
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However, from theoretical considerations it has been 
pointed out 23 that in-plane interactions are important 
for the MIT in V2O3 and and should not be treated as 
small perturbations. This view is supported by our mea- 
surements which suggest that in-plane interactions are a 
decisive ingredient to the MIT by triggering the modifica- 
tion of the V0-V1 distance via changes in hybridization. 
In this context it is interesting to note that the Debye- 
Waller factor of the in-plane vanadium neighbors does 
not display substantial changes at the transition. If struc- 
tural degrees of freedom were the driving force behind the 
increase of the in-plane interatomic distances, the Debye- 
Wallcr factor should be strongly affected. A possible ex- 
planation could be that the change in hybridization is 
driven by orbital interactions. The change in hybridiza- 
tion will also result in a change of the crystal field split- 
ting and accordingly alter the orbital occupation. Evi- 
dence for orbital interactions is given from neutron scat- 
tering results by Bao et al2&> As shown by Laad et at, 
the variation of orbital occupation can result in a first 
order phase transition^ 

Conclusion 

In conclusion, the measurements of the local structure of 
pure V2O3 in the vicinity of the metal to insulator tran- 
sition show that complex changes of the interatomic dis- 
tances on a local scale are decisive for the changes of the 
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